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Abstract. 

We study in this paper three different theories of gravitation with massive gravitons 
- the modified Fierz-Pauli model, Massive Gravity and the bimetric theory proposed 
by Visser - in linear perturbation theory around a Minkowski and a flat Friedmann- 
Robertson- Walker background. For the transverse-traceless tensor perturbations we 
show that the three theories give rise to the same dynamical equations, to the same form 
of the tensor Sachs- Wolfe effect, and consequently to the same form of the Boltzmann 
equations for the radiative transfer in General Relativity. 

We then analyze vector perturbations in these theories and show that they do not 
give the same results as in the previous case. We first show that vector perturbations 
in Massive Gravity present the same form as found in General Relativity, whereas in 
the modified Fierz-Pauli theory the vector gravitational-wave polarization modes (^'3 
amplitudes in the Newman-Penrose formalism) do not decay too fast as it happens 
in the former case. Rather, we show that such ^'3 polarization modes give rise to 
an unusual vector Sachs- Wolfe effect, leaving a signature in the quadrupole form 
Y2,±i{(^, f) on the Cosmic Microwave Background Radiation polarization. We then 
derive the details for the Thomson scattering of CMB photons for these ^3 modes, 
and then construct the correspondent Boltzmann equations. Based upon these results 
we then qualitatively show that ^'s-mode vector signatures - if they do exist - could 
clearly be distinguished on the CMB polarization from the usual ^'4 tensor modes. 

We also estimate that the graviton mass limit for the vector modes is m = 10"^^^ ~ 
10~^^cm~^ , so that vector modes with masses below this limit exhibit the same 
dynamical evolution as the massless gravitons. 

We argue at the end of this paper that CMB polarization experiments can be 
decisive to test alternative theories of gravitation by measuring CMB polarization in 
the i?-mode. 
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1. Introduction 

Among the most important basic predictions of the Theory of General Relativity (GR), 
we know that only one of them has not been yet directly tested: the existence of 
Gravitational Waves (GWs). The GWs in the GR present two basic features: first, the 
particle associated with the wave, the graviton, is massless; second, they have only two 
polarization states. The first feature comes exactly from the weak-field approximation 
in the space-time metric, yielding a spin-two massless particle field equation; the second 
feature comes from the gauge invariance of the field equations under local coordinate 
transformations. These two features are strictly related, since a massless theory allows 
gauge transformations, and such transformations lead to two polarization states. 

However, a metric theory of gravitation, as we shall see, allows up to six polarization 
modes in general [1]; furthermore, it was shown in [2] that spherical mass-resonant GW 
detectors might detect up to six polarization modes. These two facts, one theoretical 
and the other experimental, motivate us to investigate more carefully a general metric 
theory of gravitation in order to make it possible to predict further features that in the 
usual GR might, in principle, be lacking. 

An immediate generalization of GR may be constructed taking into account some 
analogies with the basic ideas of Quantum Field Theory (QFT). The simplest models 
in QFT usually involve free real massless fields with global internal symmetries as, for 
example, the electromagnetic field; then, to these models we can add further elements 
like mass, non-abelian internal symmetries, and so forth. Taking into account again 
the example of the electromagnetic field, the next generalization of this model is the 
introduction of a mass (via the Higgs mechanism), giving rise to a massive spin-one 
particle. These two vector models, the massless and the massive ones, are not theoretical 
conjectures; they give rise to observable particles, the photon and the gauge bosons 
and respectively. Then, since in GR the metric tensor field ga/s gives rise only to 
massless particles, the next extension consists in the introduction of massive degrees of 
freedom into the metric - the so-called ^''massive gravitons". There are some different 
effort towards a theory of gravitation with massive gravitons; in this paper we shall 
consider the approaches worked out in the references [3], [1], [5], [6] and [7]. 

The approach followed in references p] and [3] is an improvement of an earlier model 
pushed forward by M. Fierz and W. Pauli in the thirties [8]. It is constructed by adding 
a mass term to the linearized Einstein-Hilbert action through the prescriptions of field 
theory, that is, a quadratic term in the fields hai^ appearing in the weak-field limit. We 
shall call this model as "modified Fierz-Pauli" henceforth. In this model, GWs have six 
polarization modes [9] , which may introduce interesting features in cosmology as we are 
going to show in this paper. 

The second approach, devised by Matt Visser, [5] is based upon a bimetric theory, 
first pushed forward by the work of Rosen [10] . In this model the mass term is introduced 
by a quadratic term depending not only upon a dynamical metric as in GR, but upon a 
nondynamical background metric as well. De Paula et. al. [11] showed that this theory 
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leads to exactly the same dynamical equations as in the modified Fierz-Pauli model in 
the weak- field approximation, so that they are absolutely equivalent in this limit. They 
also have shown that GWs in this model also possess six polarization modes as expected 

In the theory of Massive Gravity, as introduced in [6], the Lorentz-invariance of 
the mass lagrangian is broken in the following way: the quadratic term in the metric 
perturbations is split into components, giving rise to five possible combinations, with 
each combination having a different coefficient; then, these five different coefficients 
are interpreted as five mass parameters, each of them being proportional to a common 
scale m. It was also shown in [6J that the Fierz-Pauli model is promptly recovered 
through a suitable choice of the mass parameters. In [7], the masses for the gravitons 
are introduced through a very clever analogy with the Higgs mechanism in QFT: the 
Lorentz invariance is spontaneously broken by a convenient choice of the "vacuum" for 
the Goldstone fields, which leads to the model devised in [6]. 

In all the models introduced above the van Dam-Veltmann-Zakharov (vDVZ) 
discontinuity [12], [13] is absent [1], [5], [6], [7|. In particular, it can be shown that 
Massive Gravity is absolutely free of ghosts and classical instabilities provided the mass 
parameters obey some constraints [7]. 

In addition to the efforts towards a direct observation of GWs (see, e.g., ref. |14j). 
we might consider, as an excellent alternative for the time being (since the GW detectors 
have not yet reached the proper sensitivity to make such direct observations), an 
indirect approach, using the recent work performed on the theory of Cosmic Microwave 
Background Radiation (CMB) anisotropics and polarization. The CMB spectrum might 
hide some "tracks" of GWs, or signatures, left by the interaction of the primordial 
cosmological gravitons with the CMB photons. The next generation of CMB satellites 
after Planck is expected to measure the CMB polarization to a better degree of accuracy, 
and it can hopefully shed a light on the investigations of primordial GWs through the 
measurements of the so-called E and B-polarization modes [15]. This means that CMB 
polarization measurements might not only be decisive to detect GWs, but it might also 
shed a light on the nature of the gravity itself; in other words, CMB measurements 
could be decisive to test alternative theories of gravitation. This is exactly our goal 
in the present work: the analysis of CMB polarization induced by GWs in a theory of 
gravitation with massive gravitons. 

To this end, the present paper is organized as follows: in Section [2] we review the 
classification of plane GWs in the Newman- Penrose formalism for an arbitrary theory 
of gravitation, deducing their six possible polarization modes. In Section [3] we review 
the three different approaches to include massive gravitons into GR, and also discuss 
cosmological perturbations in these theories. Then, in Section H] we review the basics 
of radiative transfer in the presence of weak gravitational fields, laying the ground for 
discussing the Sachs- Wolfe effect in a theory of gravitation with massive gravitons in 
Section [5l and the effect of GW vector longitudinal modes on the Thomson scattering 
in Section [61 In Section [7] we get together all the results, obtaining the related radiation 
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transfer equations (Boltzmann equations) for each polarization mode. At the end of 
this paper we discuss the obtained results and make the correspondent conclusions. 



2. Polarization States for an arbitrary Metric Theory of Gravitation 

The polarization states for a GW in an arbitrary metric theory of gravitation are given 
by the independent modes of the Riemann tensor. In order to compute its components 
in an Lorentz-invariant scheme, it is convenient to introduce, following the pioneering 
work of Newman and Penrose, [16], the quasiorthonormal complex-null basis {k, I, m, m), 
where k and I are real null- vectors and m and m are a pair of complex numbers, satisfying 
the following orthogonality relations: 

k ■ I = 0, m ■ rh = —1, k-rh = k- m = l- fh = l- m = 0. (1) 

We follow [1] and choose the following set of null vectors, 

= -^(1,0,0,1), / = --^(1,0,0,-1), (2) 

m = - -^(0, 1, z, 0), m = --^(0, 1, -z, 0), (3) 

all satisfying ([1]). With the basis given by ([2]) - ([3]) we can split the Riemann tensor 
into its irreducible parts, namely, the Weyl tensor, whose ten independent components 
are given by five complex scalars (\E'o, ^E'l, ^^2, ^3, ^4), the Ricci tensor, whose nine 
independent components are given by the scalars $00, $01, $02, $10, $20, $11, $12, 
$21, $22, and the Ricci scalar A. Throughout this section, we consider only plane GW 
propagating in the z direction, whose time dependence is given by cosut, and we use 
for the fiat Minkowski metric rjap =diag{+, — , — , — }. 

In this context, we may prove that the differential and algebraic properties of the 
Riemann tensor reduce the number of independent components to six fl], given by 

i) The Weyl tensor: 

(4) 
(5) 

(6) 
(7) 

a) The Ricci tensor: 



$02 = $20 = 0, (8) 

(9) 
(10) 
fll) 



^0 = 


^1 = 0, 


^2 = 


— ■T:Rlklk, 
D 


^3 = 


- -Rikm 


\I'4 = 


Rlfhlfhi 



$00 


= ^oi = $10 = 


$22 


Rlmlrhi 




3^ 


$11 


= 2^- 




$12 


= $21 = ^3; 
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a) The Ricci scalar: 



A = — ^'2- 
2 



(12) 



We can reduce, therefore, the number of independent components of the Riemann 
tensor to the set 

{^2, ^3, ^3,^4, ^4, $22}. (13) 

Henceforth, we call f[T^ Newman-Penrose (NP) amplitudes. They play the role 
of definite helicity states s = (0, ±1, ±2) under rotations around the z axis in a 
nearly Lorentz coordinate frame. In particular, the two real NP amplitudes (^^2, ^22) 
correspond to the state s = (which defines the scalar modes), whereas the complex 
NP amplitudes (\E'3, ^3) correspond to s = ±1 (vector modes), and {"^4, ^4) to s = ±2 
(tensor modes). These polarization modes can be represented on the x — y, y — z or 
X — z plane as can be seen through the effects of a GW on a ring of dust particles (see 
reference [TT] for further details). Now, following [Ij, it is very useful to introduce the 
"driving-force matrix" 

S,,{t) := R^ojo{u), (14) 

where t is the proper time and u = t — z/c represents a null "retarded time" as measured 
by an ideal detector in the coordinate system {t, x*}. From (fl^ we define a basis for 
the GW polarizations as follows: first, we represent the NP amplitudes as 

Pi{z, t) = ^2iu), P2iz, t) = Re ^3(m), P3{z, t) = Im ^3(m), (15) 

Pi{z, t) = Re ^^4(m), P5(^, t) = Im ^^{u), p^iz, t) = $22(m); (16) 

(from now on, whenever the polarization index r appears, it will always indicate a given 
NP amplitude according to the sequence given in fll5tll6p . so that r = 1 stands for 
"^2, and so forth). Now, writing the NP amplitudes ([5]), ([6]), ([7]) and ([9]) in Cartesian 
coordinates (recall that the NP basis given by ([2]) - (I3l) can be written in terms of the 
coordinates {t, x*}), we get the following result 

-2p2 \ 

2P3 



S 



I -\{pi+p&) 
-2p2 



\ 

6 

Pr {Z, t) Er 



(17) 



-Qpi J 



(18) 



where Ej. (z) are the basis polarization matrices, given by 
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0\ 




(19) 



Therefore, the polarization of a GW in an arbitrary metric theory of gravity can be 
fully described by the basis polarization matrices E^. {z). However, due to the tensorial 
character of the space-time metric it is convenient to cast the polarization basis (1191) 
into a tensor; hence, along with its spatial components, given by {Er)^^ (z), there are the 
00 and Oi components, which are zero by the very definition of the "full driving-force 
matrix" S ([HD, so that 



'S'oo(^) = -Roooo(^) = 0, S^iit) = RoiQo{u) 
Hence, the polarization tensor assumes the form 
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3. The Theory of Gravitation with Massive Gravitons 



3.1. The modified Fierz-Pauli model 

Let us now analyze how do GWs arise in the case of the Fierz-Pauli modified model. In 
this case, the graviton mass lagrangian appears as a quadratic term in the perturbation 
of the metric tensor hap in the weak-field limit, so that its action is given by [5], [1]: 



S 



d X 



(22) 



- AMpfh^pT'-f - ( h^ph^^ - rrh^ 
where h is given by 

h = Vaph''", 

and Mpi is the Planck mass. If instead of the contribution m^/i^/2 to the last term on 
the right-hand side of (1221) one had m?h'^ , this model would correspond to the original 
Fierz-Pauli action, which is plagued by the vDVZ discontinuity [1]. 



(23) 



CMB Polarization in Theories of Gravitation with Massive Gravitons 



7 



The Einstein equations associated with the action (1221) are given by 



- Vapd^d^h + w? {h^p - ^r]^phj = -2MpfT^p- (24) 

then, imposing the conservation of the stress energy-momentum tensor, Vq,T"^ = 0, we 
get the following constraint to the field hap on a Minkowski background, 

dji''^ = 0, (25) 

where we have defined 

hai3 = hap - ^r/o/j/i. (26) 

The equation fl2^ is exactly the same found in GR, but in the present case it 
emerges as a constraint from the conservation of the stress energy-momentum tensor 
rather than a gauge choice, as in GR. This constraint eliminates four degrees of freedom 
out of the ten independent components of the space-time metric, leaving then only six 
independent modes. Since these modes correspond exactly to the polarization states of 
the GW, we may readily associate the components of h^p with the correspondent ones 
of (!2T|) . so that the only nonzero contributions are the spatial components hij. 

Using the arguments above and plugging equation (^3]) into we obtain, in the 
absence of sources, 

(d'^d^ + m^) hij = 0, (27) 

which is clearly a Klein-Gordon equation for a wave propagating in the direction k = z. 
For the sake of simplicity we henceforth drop the bar over the tensor on the left-hand 
side of (l26l) and simply write it as hij. 

Due to the oscillatory character of equation (|27|) we may expand the tensor field 
hij into the Fourier modes as follows, 

/oo ^3 h, _ 
-oo (2^^^^' ^^^^ 

which enables us to write down the following decomposition in terms of the polarization 
tensor (EU): 



h^,{k) = Y.<,ikWik). (29) 

r=l 

In particular, for the transverse-traceless (TT) component of the tensor perturbation to 
the metric hij (corresponding to the mode with r = 4, 5), we write 

hij = et^{k)h\k)+e',^{k)h'{k), (30) 

so that Fourier transforming (!30!) to get h^j in the configuration space, we see that this 
tensor satisfies (1271) . that is 

{d^d^ + m^)hij = 0, (31) 
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which reduces to the GW equation for GR in the hmit m = 0. The tensor hjj 
encompasses then both transverse polarization modes "+" and "x" characteristic of 
GR. 

We write the extension of definition (1301) to the \t'3 modes (associated with r = 2, 3) 

as 

H, = el{k)h'{k)+el{k)h;'{k), (32) 

which corresponds to a longitudinal polarization state. As we shall see in Section |5l 
these modes can also induce an unusual angular pattern for the CMB photons and, in 
consequence, induce a different signature in its polarization pattern. 

In this paper we do not consider the GW scalar polarization modes and $22? since 
they couple to the 6goo scalar component in the metric perturbation on cosmological 
scales, and then do not produce "handedness" to excite the CMB B-polarization mode. 



3.2. The bimetric model 

The bimetric model proposed by Visser in |5] combines into a single theory a dynamical 
space-time metric with a nondynamical metric which allows the introduction of a mass 
for the graviton even in the strong- field limit. However, as was shown by De Paula et. al. 
[TT] . in the weak-field limit the bimetric model is absolutely equivalent to the modified 
Fierz-Pauli model as discussed above. Due to this property we henceforth consider the 
modified Fierz-Pauli solely. 



3.3. Massive Gravity 

In the case of Massive Gravity, as we have pointed out in the Introduction, the underlying 
ideas are a bit different. Let us now write down the space-time metric in the weak- field 
limit as 

gaf3 = Va/B + Sgaf3, (33) 

where Sgap plays the role of a perturbation to the Minkowski background metric, and 
assumes the form [6], 

^^foo = 2^9, 6goi = Si - diB, 

59^, = -Xr3 - ^^F, - d,F, + 2 (^5,, - d,d,E) , (34) 

where ip,ip,B,E are scalar fields, Fj and Si are vector fields, and Xij is a tensor field. 
The vector and tensor fields of (!34l) satisfy the well known constraints 

X./'=0, x\ = 0, F\i = S\, = (35) 

necessary to match the number of independent fields to the ten independent components 
of the metric 6 gap- Now, we construct the mass contribution to the action of the theory 
as a quadratic term in the tensor fields as well, but breaking the Lorenz invariance as 
suggested in [6] 

Mi 



''PI 



^o^9oo + 2r«i55'oi - + rn^^ga^gjj - 2m^6goo5gii . (36) 
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Now, plugging the decomposition into the mass lagrangian and in the usual 
Einstein-Hilbert lagrangian of GR, and adding up these terms, we get the following 
equation for the tensor field Xij [H] 

{d^-d, + ml) x^J = 0, (37) 

which is exactly the same equation for the TT metric contribution of the modified 
Fierz-Pauli model given by (13T]) . representing the "genuine" GW polarization modes 
for Massive Gravity. As we have argued in Section 13.11 the content of the modified 
Fierz-Pauli model plays the role of extra GW longitudinal polarization modes; however, 
in Massive Gravity, the vector fields Fj and Si evolve as massive spin-one particles with 
transverse polarization, and have nothing to do with extra GW polarization states. 

Besides this original approach devised by V. Rubakov [6], there is another way of 
breaking the Lorentz symmetry of the massive lagrangian which resembles the Higgs 
mechanism in QFT. In this case, we introduce a set of four Goldstone fields (j)^{x) and 
0*(x), i = 1,2, 3, and a scalar, a vector and a tensor field constructed as follows [7j, [T8] : 

X = A-^(?"^a„0°a^0°, (38) 

r = A-^^7"^9,/9;30\ (39) 

W'^ = A-^g'"^d^(f)'df^(f)^ - X-^V'V^, (40) 

where A is the energy scale of the theory. With such elements we introduce an arbitrary 
function F = F {X, V^, W^^), so that the action for Massive Gravity reads 



S 



J d^xV^ [-Mj,iR + A^FiX, V\ + Cmatter] , (41) 

where the first term on the right-hand side represents the usual Einstein-Hilbert action, 
and Cmatter IS the lagrangian for ordinary matter minimally coupled to the metric. We 
can show that in the linear perturbation regime given by (l33l) . after setting the Goldstone 
fields to their "vacuum" values, 

9af3 = Val3, <f = AH, 0^ = AV, (42) 

and substituting these values and the decomposition flM|) into (HTi) we get exactly the 
same mass lagrangian as fl36l) . with the mass parameters now being related to the 
functions F and their derivatives [7j, [18]. 

3.4- Cosmological Perturbations for Massive Gravitons 

Once we have discussed the key features of metric perturbations around a flat Minkowski 
background, let us now address the same question in a flat Friedmann- Robertson- Walker 
(FRW) background. We start with GR in which, in the standard theory of cosmological 
perturbations [19], the metric 6ga/3 is decomposed exactly in the same way as we did in 
but, in this case, we multiply all these components by the square of the scale factor 
of the universe 0(77) (r/ is the conformal time), that is 
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The constraints for the cosmological case are the same as (15^ . Now, since we are 
interested in the CMB polarization induced by GWs, we focus only on the TT part of 
the metric perturbation; in this case, the Einstein equation for the tensor field is given 
by (the prime ' indicates a derivative with respect to conformal time) 

h",, - V^hij + 2nKj = 0, (44) 

which simply describes transverse GW travelling on an expanding background. 

In the case of the modified Fierz-Pauli model the same metric decomposition cannot 
be performed due to the extra polarization modes; we instead introduce 

where (j) and Q are scalar fields, Xi is a divergenceless vector field, and hij is the 
cosmological version of the tensor given by the solution to equation fl271) . carrying the 
correspondent six polarization modes spanned in the NP formalism. The two scalar 
fields, plus the two components of the transverse vector field and the six modes of the 
tensor field give exactly the required ten degrees of freedom. The mass lagrangian for this 
model can be constructed analogously as in fl22l) . that is, it appears as a quadratic term 
in the metric fHSl) . The full action is then obtained by adding up this contribution to the 
usual Einstein-Hilbert one, and the Einstein equations can be derived using the standard 
tools. Before doing that, it is convenient to decompose the tensor perturbation hij into 
its TT and longitudinal parts in the Fourier space. In (1281) the whole time-dependence 
of hij is contained in the exponential since it is a solution to a wave equation of the 
form ( !27ll ; now, such time-dependence changes because of the extra temporal function 
a{r]) appearing in fH5l) . which introduces a damping in the oscillation. Therefore, we 
Fourier-expand the massive tensor perturbation hij as 

= / 77r^hi,{r^,k)e-^'^-\ (46) 



where 



-oo (27r)3/2 



kAv,^) = E4(k)^'(^,k), (47) 



so that the TT and longitudinal components of hij can be written in the same foot as 
(I30D and ([32]), that is 

hij (r/, k) = 4(k)/.^ (r/, k) + 4(k)/.5 (r/, k) , (48) 

4 iv, k) = elik)h' (r^, k) + 4 (k) h' iv, k) ; (49) 

now, extracting the Einstein equations from the action for the cosmological Fierz-Pauli 
model as we have sketched above, we see that both fields hjj and hjj satisfy the same 
dynamical equations 

h^'\j - V^h^.j + 2nhij + a^m^hij = 0, (50) 

h\\" ij -V'^h^ij + 2n}iij + a^m^h}lj =0. (51) 
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In the case of Massive Gravity, the cosmological perturbations to the metric are 
given by fH3|) . together the following set of perturbations to the Goldstone fields in the 
unitary gauge (H2l) . 

= 00 ^ ^2^0^ 0^ = 0^ + (a^ + A'*) , (52) 

where A" e A are scalar fields and A* is a divergenceless vector field. Now, under 
infinitesimal coordinate transformations 

f] = V + e, S:' = x' + e, (53) 
we can show that the following vector fields 

Wi = Si + F[, cTi = Xi- Fi, (54) 

are invariant. 

The action for Massive Gravity on a fiat FRW background is then given by (l4Tll 
with fH3|) and the Goldstone fields set to their vacuum values fH2l) ; the matter lagrangian 
'^matter IS assumed to be described by a perfect fluid whose perturbations for the fluid 
four-velocity are 

6ui = a{Ci + di(), 6uq = aip. (55) 
With these features, the Einstein equations for the tensor field Xij are given by [20] , 
X'\j - V\., + 2nx[, + a'^mlx^, = 0, (56) 
whereas for the gauge-invariant vector fields defined by fl54l) the Einstein equations read 
\/^zUi-2a'^PmMpi^{l + w)Ci = 0, zu- + 2nzui- a'^mlai = 0, (57) 
m^VVi = 0, (58) 

where St^ = ( — {E' + B), and w is the parameter appearing the equation of state of the 
ordinary matter, p = wp [20] . 

Solving equations (1571) - (l58ll we conclude that the only relevant vector field is Wi, 
whose amplitude decays with [20j, which is exactly the same behavior of vector 
fields as derived in GR. 

Therefore, the Fierz-Pauli modified model and Massive Gravity give rise to the 
same results for the TT polarization modes of the tensor perturbations as can be seen 
from equations (|50l) and (1561) . whereas for vector perturbations the situation changes 
drastically. In the modified Fierz-Pauli model the vector modes of GW polarization 
obey the same equation as the TT modes, ( l5Ti) . so that they really may contribute to 
the polarization of CMB as we have shown in [21J for the field Xij'-, however, in Massive 
Gravity, the vector perturbations behave exactly as in GR, which means that they decay 
too fast after the inflationary phase and do not leave any signature on CMB polarization. 

Therefore, the vector modes of the modified Fierz-Pauli model, unlike the 
predictions of GR, Massive Gravity, and other modified models of gravity, can leave 
signatures on the CMB. To see how this can be achieved, let us estimate now a limit 
mass to be detected by CMB polarization experiments, according to the discussion in 
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pi] . Since the equations for the tensor (I3U]) . and vector modes are identical, the 
dispersion relations are also identical, given by 

uj^ = k'^ + m^; (59) 

now, since in GR only GW with frequencies u within the range 10^^^ Hz to 10~^^ Hz may 
leave a signature on CMB polarization [22|, we see that these frequencies correspond to 
comoving wavenumbers k within the range 10~^^cm~^ to lO^^^cm^^ . We then use the 
values of k of GR, varying the frequencies in order to obtain constant nonzero graviton 
masses through equation ( l59i) . As a result, if the values of the graviton mass m lie within 
the range 10~^^ - 10"^^^', the correspondent frequencies have values very close to the 
expected in GR; in particular, there is a graviton mass limit, m = 10~^^g ~ 10~^^cm~^, 
so that below this limit the dynamical evolution of the massive modes is indistinguishable 
from the tensor massless modes. Hence, since the vector modes in the modified Fierz- 
Pauli model are governed by the same equation as the tensor modes, then the graviton 
mass limit for the vector modes should be exactly the same. 

Along with the value discussed above for the vector massive modes, there are other 
mass limits obtained so far in the literature by using different tests. For instance, 
Goldhaber and Nieto [23] found a limit m < 2.0 x 10~^^g analyzing the motion of galaxies 
in clusters. Later on, Talmadge et al. [M] studied the variations of Kepler's third law 
when compared with the orbits of Earth and Mars, and found a limit m < 7.68 x 10~^^(y'. 
Recently, Finn and Sutton |3j calculated the decay of the orbital period of the binary 
pulsars PSR B1913+16 (Hulse and Taylor pulsar) and PSR B1534+12 due to emission 
of massive gravitons, and found m < 1.4 x 10~^^g. 



4. The Radiative Transfer Equation in the presence of Weak Gravitational 
Fields - an overview 



Once we have discussed the key features of the different theories of gravitation 
with massive gravitons, we now turn our attention to the polarization of CMB. To 
do so, we initially review the theory of the radiative transfer in the presence of 
weak gravitational fields, following closely the seminal paper by Polnarev [25j (for a 
pedagogical introduction, see [26]). 

Let us consider a given beam of radiation characterized by its Stokes parameters 
{Ii, Ir, U}, where /; and Ir are the intensities of the radiation in the directions / and r, 
respectively; / = + is the total intensity of the wave, and the parameter Q is given 
hj Q = Ii — Ir- These functions are strictly related to the photon distribution function, 
which can be cast in a symbolic vector of the form [27] 



f 



Ir 



(60) 



where U is the other parameter associated with linear polarization. The components of 
/ are functions of the conformal time rj, the comoving spatial coordinates r, and also of 
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the photon angular distribution. Since the photons are scattered by the free electrons 
prior to recombination via Thomson scattering, their distribution function will be shifted 
according to the equation [25] 



df ^p" df ^ df dv _ ^ - 
drj p dx'- dv drj ' 



(61) 



where v is the photon frequency, p is the photon momentum, C[f] is the scattering term 
given by 

1 n 



C[f] = -arN^a 



f 



d^i'dif P {fi, ip, fi', v?') / 



(62) 



(Jt is the Thomson scattering cross-section, is the number of free electrons in the 
unit comoving volume, /i = cos 9, and P {fi, ip, fi', ip') is the scattering matrix given by 



P(/i,<^,/i',<^') =P 



p' (/i, /.') + Ji- fiyi - fi'^p' (/., ^, /i', ^0 



(63) 



where 



P 



/ 1 \ 

10 

10 

V 2 / 





( 2(1 


-/.2)(i-/i'2) + /.V2 







\ 


3 




/^'^ 


1 








4 


















V 











flfi' J 



p' 



3 
4 



f ififi' cos ip 



2 fj,' sin ip 

\ 



-2/i siiaip \ 



cosip 

cosip J 



(65) 



/ /i^/i'^ cos 2V^ — /i^cos27/' — /i^yu' sin 2?/' \ 



— /i'^ cos 2ip 
fifi''^ sin 2'ip 




cos 2ip 
— /i sin 2ip 




fi' sin 2il) 
fifi' cos 2?/' 






0/ 



(66) 



where we have defined ip := ip — ip' . 

In equation fl6T|) the third term on the left-hand side shows the influence or signature 
of the GWs on the CMB polarization, whereas the term on the right-hand side gives 
the details of the Thomson scattering. Since the term (!62l) depends upon the photon 
angular function through the angles {fi, ip), we see that different GW imprints may lead 
to different processes of polarization via Thomson scattering. In the next two sections 
we address these issues in detail. 
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5. The Sachs- Wolfe effect induced by Massive Gravitons 

As we have seen in Section 13.41 the modified Fierz-Pauli model and Massive Gravity are 
equivalent with respect to the GW TT tensor modes, but only the vector modes of the 
first model may give rise to relevant contributions to CMB polarization. Hence, from 
this section on, we shall always consider only the Fierz-Pauli model to address all the 
issues concerning GW vector modes. 

We have discussed in Section H] that the CMB photons are polarized due to the 
Thomson scattering with the free electrons in the epoch of recombination. Prior to 
Thomson scattering, the cosmological perturbations imprint a signature on the photon 
angular pattern, the so-called Sachs-Wolfe (SW) effect, which can be understood as the 
shift of photon frequency along the line of sight. 

This effect can be computed through the geodesic equation for the photon; then, 
using the metric perturbation fH5l) around a fiat FRW space, and the constraint 
gai3P°'p^ = 0, where is the photon four-momentum, the geodesic equation reads 
for the tensor field hij, 



du 



where A is an afiine parameter. The product e^p'^p' appearing in fl67|) may be evaluated 
as follows: first, since the photon travels along an arbitrary direction p, let us construct 
a reference frame around the GW propagation vector k such that the GW polarization 
modes assume the simplified form given by fl21l) . To this end we introduce the 
polarization vectors {^(i),^(2)}5 defined by 

^ij = - ^(2)i^(2)i5 (68) 

where e^j is given by (pT!) for each polarization component r, satisfying 

^(1) ■ ^(2) = ^(1) ■ k = £(2) ■ k = 0. (69) 

The trihedron 

{^(i),^(2)>k} (70) 

is then a reference frame around the GW direction k in which the polarization tensor 
preserves the simple form (12T!) . 

Second, let us express the photon momentum p in terms of spherical coordinates 
around (1701) . so that the following relations hold, 

k ■ p = cos 9, • p = sin 9 cos Lp, £(2) ' P = sin 9 sin ip, (71) 

therefore, using (!2T|) and (17T|) . it follows that 



4pV = fiJl- fi'^cos^ (X F2,+i (/^, <^) , (72) 



^ijPY = - fi'^siuif (X ¥2,-1 (/i, if) , (73) 



4pV = 1 - /i' cos2y. cx r2,+2 (/i, ^) , (74) 
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4pV = (l - /i^) sin 2^ (X 1^2,-2 (/i, ^) ; (75) 

where Yim (/x, V') are the usual spherical harmonics. The results ( |72l) - (1751) show us that 
the GW imprint upon the photon angular distribution is in the form of a quadrupole, 
with m = ±2 for the \E'4 modes (r = 4, 5, which coincides with GR), and with m = ±1 
for the \E'3 modes (r = 2, 3). 

Hence, from relations P6|) . (H7I) . fl671) and (!72|) - (175|) we get the following geodesic 
equations: 



1 C/Z/Q 

« ^^2,±1 (/^, V') , (76) 



a) "$4 



Uq drj drj 
IdVo l9/i4'5 

— -1- « -77^^>^,±2 (/i, f/') (77) 
Vq drj 2 a?] 

where z/q = va{ri). 

Then, massive gravitons with the \E'4 polarization modes give rise to the usual tensor 
SW effect in GR, whereas the \&3 modes do not yield its well-known vector version in 
GR (see, for example, the discussion in |28]). This happens because \E'3 modes arise as 
GW longitudinal states of polarization, and not as a massive vector fields as in GR or 
Massive Gravity. Anyway, they will leave a different signature on the CMB polarization 
by means of Thomson scattering, as we discuss in the next section. 

6. The Basis for Thomson Scattering 

We now turn to the derivation of the Thomson scattering term for massive gravitons. 
Let us consider first that the incident radiation prior to be Thomson scattered is 
unpolarized, with the angular pattern dictated by the SW effect (1761) and ( 1771) . In 
this case, the Stokes vectors for the incident radiation are given by 
i) ^3 : 



a) 



d^ = -/i-y/ 1 — /i^ cos (f u, d^ = -/i-y/l — /i^ simp u, (7^ 



^ (l - /x^) cos 2lp u, = ^ (l - fx'^^ sin 2lp u, (79) 



d' 

2 

where we have defined 

/1\ 
u = 1 

Now, defining the operator P as 



(80) 
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it is straightforward to see, for ^ = (f 

(82) 
(83) 



where P is the scattering matrix 
(r = 2,3,4,5), that 

Pa'' = pa'' + qif, 

where p and q are constants, and If is a basis vector such that 

py- = p'd'' + q'V, 

where p' and q' are constants as well. From (ITHl) - (l83l) we readily see that 
i) : 

—fi COS ip 

^ 2 simp J 



/ yusincyj \ 
— sin ip 
—2cosp> I 



M) 



ii) 



(1 + yU^) COS 2^9 
- (1 + /i^) COs2(y9 

4yU sin 2^9 

The basis vectors given by fl78l) and ([1 



(1 + /i^)sin2v? 
- (1 + /i2)sin2v? 

— 4/i cos 2(y9 

Bl) and fl85l) allows us to factor out 



5) 



the angular dependence of the photon distribution vectors, and they constitute the 
first-order contribution to /. The remaining contribution for / comes from the h = 
solution, that is, in the absence of GWs, in which the photon distribution vector is given 
by 

/ := /o = /o(^^o) fl. (86) 
Now, the linearized photon distribution vector can be written as 

r ~ fo + e'''-' [«'■ iv, f^, ^o) + iv, f^, ^o) b^] , (87) 

where a'' {rj, /x, z/q) and {rj, /i, z/q) are functions to be determined as solutions of the 
Boltzmann's equations (I^Tj) . 



7. The Complete Boltzmann Equations 

Now, once we have obtained the form of the photon distribution vector (IHTl) for the TT 
and longitudinal GW mode, we are able to write down the full Boltzmann equations 
They are given by 
z) ^3 : 

x" + (g - ikiJ,)x'' 



+ /3''(-l-/i'2 + 2/i'^ 



^9) 



for r = 2, 3, and 
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a) \i/ 



4 



C' + {q-^kf,)C = H\ (90) 



(91) 



lb J-i 

for r = 4,5. The functions and x ^"^^ defined as 

f = a" + /3^ (92) 
X^ = a^-^\ (93) 

and the scattering rate q is defined as g = axN^a, and finally 

H' = 1^. (94) 

The equations (pOl) and (!9T!) for massive gravitons with \&4 mode are identical to 
the correspondent ones in GR; the difference here lies on the Boltzmann equations for 
the modes, fl88|) and fl89|l . which do not appear in GR. Since the vector and tensor 
modes satisfy the same dynamical equation, and the mathematical form of the equations 
( I88|) and (189!) is very different from (pOl) and (19T1) . it is clear that the vector polarization 
modes of massive gravitons leave a characteristic signature distinguishable from the 
tensor one, which could, in principle, be probed by measurements on the CMB E and 
i?-modes. Since the experiments in the Planck satellite will improve the WMAP5 results 
for the E-mode, we may expect that such future measurements might decide whether 
nontrivial GW signatures - as we showed here through equations (!88l) and (!89|) for \E'3- 
modes - appear or not in the CMB polarization spectrum. In this case, we conclude 
that CMB polarization measurements may be decisive to test alternative theories of 
gravitation - in particular, the massive model as we discussed here. 



8. Conclusions 



We have analyzed in this work three theories of gravitation with massive gravitons, 
and we have shown that the modified Fierz-Pauli model coincides with Massive Gravity 
for the TT tensor field in linear perturbation theory around Minkowski and fiat FRW 
backgrounds. We have shown that such tensor perturbations associated with massive 
gravitons give rise to the usual tensor SW effect, lead to the same Boltzmann equations 
for these modes in GR. 

We have deduced the dynamical equations for the GW vector longitudinal 
polarization modes (^^a-modes) in the modified Fierz-Pauli model and shown that they 
do not give the same results of Massive Gravity, in which vector perturbations behave like 
in GR; instead, they give rise, in a cosmological scenario, to a nontrivial SW effect which 
leaves a vector signature of the quadrupolar form Y2^±i{fi, (p) on the CMB polarization. 
Also, such massive vector modes possess a mass limit of m ~ 10~^^cm~^ as the tensor 
modes, so that below this limit the Fierz-Pauli model is absolutely indistinguishable 
from GR in terms of the dynamical evolution. 
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Analyzing the Einstein equations for such \l/3-modes we concluded that these vector 
signatures could be present at recombination epoch, unlike the vector perturbations 
in GR and Massive Gravity, which would decay too fast and would not leave any 
signature on CMB polarization. Therefore, we calculated the new basis for the Thomson 
scattering for such ^3-modes, and then deduced the appropriate equation for the 
radiative transport. Based upon these results we have shown quahtatively that ^3-mode 
vector signatures could clearly be distinguished on the CMB polarization from the usual 
tensor modes if the former do exist; hence, we could look for such signatures in the 
£'-mode performed by Planck satellite. 

In this sense we argued that Planck polarization measurements could be decisive 
to test alternative theories of gravitation. 
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